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Abstract 

A simple banking network model is proposed which features multiple 
waves of bank defaults and is analytically solvable in the limiting case of an 
infinitely large homogeneous network. The model is a collection of nodes 
P^ , representing individual banks; associated with each node is a balance sheet 

f ri ■ consisting of assets and liabilities. Initial node failures are triggered by 

external correlated shocks applied to the asset sides of the balance sheets. 

These defaults lead to further reductions in asset values of all nodes which 

in turn produce additional failures, and so on. This mechanism induces 

indirect interactions between the nodes and leads to a cascade of defaults. 

^ , There are no interbank links, and therefore no direct interactions, between 

OO ' the nodes. The resulting probability distribution for the total (direct plus 

r*~ I systemic) network loss can be viewed as a modification of the well-known 

[^ . Vasicek distribution. 



The purpose of this short note is to introduce a banking network model 
capable of describing cascading bank failures. Consider a banking system where 

^o I each bank is represented as a node. Each node i has an associated balance sheet 

which consists of assets Ai and liabilities Li. Nodes are solvent if Ai > Li, and 
they default when this inequality no longer holds as a result of the stressed asset 
side. All nodes are assumed to be initially solvent. The first wave of defaults is 

/\ . triggered by external correlated shocks applied to the asset sides of the balance 

sheets. The shocks change the asset values from Ai to new values Ai^i. Some 
of them go below the liabilities {Ai^i < Li), and the corresponding nodes fail. 
These defaults lead to extra shocks for assets, changing the asset values from 
Ai^i to Ai^2- As a result, some of the nodes which survived the initial shocks 
fail to satisfy the new solvency condition Ai^2 > Li, and default. The second 
default wave produces yet another set of shocks for asset prices, more nodes fail, 
and so on. All but the initial shocks are modelled by discounting asset prices 
Ai^i, and the discount factor applied to Ai^i after k default waves has the form 

I?fc =exp(-agfc) (1) 
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where qk is the network loss at this stage defined as the fraction of defaulted 
nodes, and a is a constant [T][3]. Asset prices Ati are assumed to be lognormal; 
the asset log changes (returns) 

i?.^ln^ (2) 

are jointly normally distributed with identical means /i and covariances <j'^ p, 
where a is the standard deviation and p is the correlation coefficient. These 
variables admit the representation 



Ri ^ p + (TXi, Xi = y^Z + ^1 - psi (3) 

where Z is the market factor common to all nodes, and Ei are the node-specific 
idiosyncratic factors. Variables Xi are jointly standard normal with correlations 
p, while Z and all Si are standard normal and mutually independent (e.g., 0]). 
It is convenient for what follows to rewrite Eq. ^ as 

Ri ^a£i + P (4) 

where 

a = a^Jl — p, (3 — p + (j^/pZ (5) 

The rest of the exposition is focused on the limiting case of an infinitely large 
(n — >■ oo, where n is the number of nodes) and homogeneous {Ai = A and 
Li — L for any i) network in which the full analytical solution is possible. 

The equation describing the cascade of node defaults in this limiting case 
can be derived as follows. The node i fails in the initial default wave if Ai^i < L. 
This corresponds, via Eqs. Q-®, to 

e^<^l, 6i^^ (in J -13] (6) 

The network loss qi at this stage, which is the direct loss caused by external 
shocks and contains no systemic contribution, is equal to the probability of the 
outcome Ai^i < L. Since Si is a standard normal variable, the fraction qi is 
given by 

q^ = P {A,^i < L) ^ P {e, < 5i) ^ N {5i) (7) 

where P {E) denotes the probability of the event E, and N (x) is the standard 
normal cumulative density function (CDF). The idiosyncratic loss q is equal to 
the probability of an individual default and corresponds to vanishing correlation 
p = 0. The relation between the direct loss qi and the idiosyncratic loss q is 
given by 

., . '^-'^p.f (8, 

together with Eq. ([T]). The first wave of defaults changes asset values from Ai^i 
to 

A,^2 = DiA.^i = Ai^i exp {-aqi) (9) 



triggering the second default wave. The network loss increases from qi to 

q2 = P{A^2<L) = P{e,<S2) = N{S2) (10) 

where 

S2 = -(lnj-/3 + aqi] =Si + KNiSi) (11) 



and 

a a 



a (7a/1 — p 
After the second default wave the asset values become 



(12) 



A,,3 = Da^,,! = Ai^i exp {-aq2) (13) 

and the network loss increases to 

q3 = P {A,3 <L)=P{e^< (5,3) = N {6z) (14) 

where now 

<53 = Vln ^ - /? + a(?2 j =5i + nN (Ja) (15) 

This process continues to infinity. After k default waves the network loss is 

qk^P (A,,fc <L) = P (£, <5k)=N {6k) (16) 

and the cascade equation has the form 

Sk = F{Sk^i), F{x)^5i + KNix) (17) 

Eq. (flTl) is a one-dimensional iterated map, the set {x, F {x) ,F {F {x)) , . . .} is 
called the orbit of x under F, and x is the initial value of the orbit (e.g., [HIS]). 
Function F (x) is continuous and increasing, so Eq. (fT7|) is an invertible map 
and can only have fixed points. When more than one fixed point exists, they 
are alternatively stable and unstable, and the unstable fixed points are the 
boundaries that separate the basins of attraction of the stable fixed points (e.g.. 

The total loss goo is the fraction of failed nodes after the cascade of defaults 
exhausts itself, i.e., after the infinite number of default waves, 

goo = N (Soc) , Soo = lim Sk (18) 

fc— >oo 

and Soo satisfies the fixed-point equation 

x = F{x) (19) 

which can be rewritten as 

f{x) = Si, f{x) = x-KN{x) (20) 



Function / (x) is monotonically increasing for k, < kq, where 

Ko = V2^ ~ 2.5066 (21) 

so in this case there exists only one fixed point for any value of di. In the 
opposite case k > kq, however, / (x) has a minimum at xq and a maximum at 
xi = —xq, where 

a;o = J21n— (22) 

V 'So 

Consequently, in the interval 

yo<Si<yi, yo^fixo), Vi^fixi) (23) 

there are three fixed points, zi, Z2 and Z3, such that 

zi < xi < Z2 < xq < Z3 (24) 

Fixed points zi and Z3 are stable, while Z2 is unstable, since 

F'(zi)<l, F'(Z2)>1, F'{zs)<l (25) 

(the prime symbol denotes the first derivative of a function), so the choice is 
between zi and Z3. The basins of attraction are (— cxd, Z2) for zi and (z2, +00) 
for Z3; the initial value of the orbit is Si, and Si < zi < Z2, which means that 
^1 belongs to the basin of attraction of the leftmost fixed point zi. As a result, 
Soo — zi for the above interval of Si values. The overall picture is as follows: for 
any Si < yo, there is a single fixed point; when Si reaches yo, a stable- unstable 
pair of fixed points, Z2 and Z3, is born at a; = xg (this event is called the fold 
bifurcation, e.g., 8,), but S^o — zi until the value Si ~ yi is reached. At this 
level of ^1, fixed points zi and Z2 collide and annihilate each other at x = cci 
(another fold bifurcation). Function g{x), linking the total loss q^o and the 
direct loss qi via the relation between Soo and Si, 

Soo =9 (Si) (26) 

jumps therefore from zi = Z2 = xi to the value z^ — X2 j^ xi of the rightmost 
fixed point which is found from the condition f {X2) = j/i- Functions gk (x), 
defined as 

Sk - gk (Si) (27) 

are, on the other hand, continuous, increasing and can be inverted, leading to 

Si = hkiSk), hk{x)=gj:^{x) (28) 

Functions hk (x) are also continuous and increasing; together with their limiting 
function 

h (x) = lim hk (x) (29) 



they are used below to calculate loss distributions for q^ and qoo ■ Function h (x) 
coincides with / (x) for k < Kg and is constructed in the case k > kq by replacing 
the segment of / (x) between Xi and X2 by the horizontal line connecting points 
(xi,2/i) and (0:2,2/1) (this is the only possibility for the limit of a sequence of 
increasing functions which includes both these points). 

The loss after k default waves qk is, via Eqs. ([5]), p^ and (l?7|) . a function 
of the market factor Z which is a standard normal variable. The probability 
distribution for qk is found by deriving an inequality for Z equivalent to the in- 
equality qk < X. Since A^ {x) in Eq. P^ is an increasing function, the inequality 
for Sk reads 

dk < N-^ (x) (30) 

Function hk (x) in Eq. (j28l) is also increasing, so the inequality for 5i is 

,5i < hk {N-' (x)) (31) 

Finally, because of Eq. ([5]) , the equivalent inequality for Z is 

1 



Z>-Ak{x), Afc(x) = — y/r^hk{N-^{x))-N-^{q) (32) 

As a result, the CDF for qk is 

Fk (x) = P {qk <x)=P{Z> -Ak (x)) = N {Ak (x)) (33) 

and the corresponding probability density function (PDF) has the form 

^(Akix)) 



p,{^x)=Fl{x) = ,J^h'^{N-^{x)) 



^{N-^{x)) 



(34) 



where (x) — N' {x) is the standard normal PDF. 

The probability distribution for the total loss qoo is obtained as the limit 
of this result when fc — > cx) and corresponds to using h (x) instead of hk (x) in 
Eqs. (I5^ - (|M)) . For the purpose of completeness, 

i^oc (x) - P (qoo < x) = iV (Aoo {x)) (35) 



..(.)^^^(.)^/^/.'(^-(.))l^l^ 



(36) 



Aoo(a;) = — yn^h{N-^{x))~N-^{q) (37) 

From the properties of h (x) it follows that in the case n > kq its first derivative 
h' (x) vanishes for xi < a; < X2 and is discontinuous sX x — X2 (it is continuous 
at X = xi, since /' (xi) = 0). This translates to the PDF poo (x) which in this 
regime is split into two parts (it is equal to zero when A^(xi) < x < A^(x2)) 
and has a jump at x = A^ (2:2)- The first of these features is consistent with the 



inability of Soo to have any value in the interval (xi,a;2), which is reflected in 
the jump of g (x). 

The probability distribution for the direct loss qi, on the other hand, corre- 
sponds to hi (x) — X and has the form 

Fi (x) =P{qi<x)^N (Ai (x)) (38) 



„,.,..; ,.)^,/^^» (-) 



Aiix)^— ^/T~pN-'{x)-N-'{q) (40) 

VP '- -' 

This is the well-known Vasicek distribution 9 . Because of the similarity in 
the mathematical structure of the two distributions and the fact that Eqs. ([55)) - 
(PO)) are obtained from Eqs. ([55| - ([57)) in the limit k — >■ 0, the derived total 
loss distribution can be considered as a modification of the Vasicek distribution 
which takes into account the systemic component of the network loss in the 
proposed model. 
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